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We compute for the first time full elastic deformations, as well as length, of self-trapped electronic 
states in carbon nanotubes of general radius and chirality, within the unifying framework of a 
recently introduced two field model for electromechanics of carbon nano-structures. We find that 
deformations are highly non monotonic in the chiral angle, whereas the length of the polaron is not. 
Applications include nano-mechanical devices as electrically or optically driven nano-actuators. 
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Coupling betw een e lectronic and elastic structure of 
carbon nanotubeJ^I^E! leads to theoretically interesting 
and technologically relevant phenomena: among other 
things, gap opening can be induced via straiii^ (allow- 
ing mechanical control of the conductance) and i nje cting 
electron/holes induces mechanical deformationaSd (im- 
portant for realization of electrical nano-actuators) or 
shifts in Raman phonon modes^^^*'. Novel excitations 
leading to spontaneous distortions driven by electron- 
phonon couplings are of parti cular interest to the elec- 
tromechanics of nano devices ^i3l2\ particular light 
induced mechanical deformations, or optomechanical ef- 
fects can be observed at high temperatures in nan- 
otube, as excitons are more strongly bound^'^. Recently, 
Verissimo-Alves et have predicted the existence of 
self-trapped electron and hole states in semiconducting 
carbon nanotubes by feeding parameters from ab initio 
density functional calculations into a very simplified con- 
tinuum model. They found lengths of the order of 40-60 
nm and energies of order of 10^^ — 10^^ meV. Their long 
polaron approximation is confirmed by a atomistic nu- 
merical analysis e.g. of Bratek et al}^. 

Chiral indices are known to play dramatic role in the 
physics of carbon nanotubes. Not only they dictate the 
metallicity, they are known to control in a v ery de licate 
way electromechanical and vibrational effect^^'^'^. It is 
of outmost importance both theoretically and technologi- 
cally to be able to predict for which chirality a polaron in- 
duces e.g. contraction rather than elongation, radial ex- 
pansion rather than shrinking, or torsion et cetera. While 
Re f.'^ first suggested the presence of polarons, their sim- 
plified model cannot address the complex deformations 
in carbon nanotubes and completely ignores chirality de- 
pendence. On the other hand, more complete treatments 
based on atomistic models are too complex to be solved 
analytically or too computationally costly to return de- 
tailed answers for a general case. 

In this letter we employ a recently introducecP^ bicon- 
tinuum model to solve this predicament and calculate the 
full range of deformations induced by self-trapped elec- 
trons/holes for the general carbon nanotube. We find 
that deformations are, as one would expect, highly non 
monotonic in the chiral indices, even with abrupt change 
of signs in nanotube of the same chiral index and about 



the same radius. 

In a recent work, Nisoli et aW^ have introduced a two 
field unifying framework for elasticity, lattice dynam- 
ics and electromechanical coupling in carbon nanostruc- 
tures that accounts for the full atomistic detail of the 
graphenic lattice, and explains a wealth of experimen- 
tal and numerical results without computationally in- 
tensive atomistic treatments. They defined an elastic 
field for each of the mutually interlaced triangular sub- 
lattices that make the honeycomb lattice, u'^{x), v'^{x), 
i = 1,2 {see Fig. fl]) and the corresponding strain ten- 
sors'^ u*^ = d^^u^\ v"^^ = d^^v^\ and wrote the elastic 
energy via considerations of symmetry. Unlike Ref.'i^, we 
work here with the average displacement 2p* = (u* -I- v^), 
inner displacement 2g* = (u* — v'), and corresponding 
strain tensors p*^ , q^^ . In these new variables, the elastic 
energy of Ref.li^ becomes Wei — J Vci[p, q] da;^, with 

+ ^C.2g2_2/3e,,fc(Jy^ (1) 

where we have employed the long polaron approximation 
to neglect the dispersion of the optical branche^^l]^ ^j^g 
tensor Cijk is invariant under the C^y group and can be 
represented by the three unit vectors {e^'''};=i_3 of Fig. [I] 

e^.^='Eef^e«e«. (2) 
1=1 

There are thus 4 parameters: tor, the graphite-like op- 
tical frequency; /3 which determines the strength of the 
rotational symmetry breaking, contains all information 
about the point group symmetry of graphene, and de- 
fines the important length £ = 4/3/a;p = 0.3^ and the 
generalized Lame symbols fi, A, which can be expressed 
in terms of the longitudinal and transverse speed of sound 
in graphene, = 2/i-|- A — 4/3^/wp, = /t— 4/3^/a;p, and 
are related to the actual Lame symbols of graphene fir, 
Xr via fir ^ fl - 4/3^ /wp and Ar = A -I- 4/3^ /wp. We then 
write p in term of the isotropic (70), anisotropic (7), and 
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shear/torsional (77) strain in the nanotube coordinates 




= 7o + 7 

lo-l 
p4>z ^pz4> 



(3) 



In these new variables the elastic energy ([T]) now reads 

^el = 2A 72 + 2/i + 72 + r^^) + g2 

- 4/3 9-^(733 -7703) + 4/3 <7^(7C3 + 77S3). (4) 

(We have shortened C3 = cos(36'c), S3 = sin(30c), is the 
chiral angle of the nanotube, and used e^^^.^ — — e^^^^z = 
-S3, and ez,z,z = -e,/, = -C3.) 

As explained in RefP^, one can simply "wrap around" 
the elastic energy of graphene to deal with carbon nan- 
otubes. In the cylindrical geometry, with coordinates 
{r, (/), z} of Fig. [1] a minimal coupling between the tan- 
gential displacements and the radial appears in 
Vol of (fTl) via p'l"!' = {d^p'f' + p") /7E3. As explained. 



we assume no azimuthal dependence, and thus 27^ = 
(d^Pz+Pr/r), 27 = {-dzPz + Pr/r), 2-q = d^p^,- We 
are neglecting the breaking of the hexagonal symmetry 
brought upon by the chiral vector that defines the wrap- 
ping of the carbon nanotube. This symmetry breaking 
allows for new terms to be introduced in Vd as curva- 
ture corrections. The parameters of our elastic energy 
for graphene are also corrected by curvature. 

As for the low energy electronic excitations in a semi- 
conducting nanotube, they can be described in terms of 
an envelope wave- function i/P^ of energy density 

E,[i},p,q]^-i,^—^TP + E,^[^,p,q]. (5) 
2m 

i?op, the coupling between phonons and an injected elec- 
tron or hole, at lowest order both in the in-plane elastic 
fields and in the electron probability density ji/jp can be 
deduced via considerations of symmetry 



I^AP [f {(fie + 7C3 + '753) + i^'7o] 



(6) 



(e is the bond length). Physically, the term proportional 
to V emerges in dielectric tubes from a deformation- 
induced change in the bandgap, whereas the term propor- 
tional to v' comes from the shift in energy at the K point 
of the Brillouin due to second-nearest-neighbor atom^. 
In the context of a simple tight-binding treatment 2, hop- 
ping integrals are modulated by in-plane elastic deforma- 
tions.fi^ di^'' = — r ef^ej'^'-' -\-Tei\^/e are the three 
nearest-neighbor hopping integrals along the three bonds 
of unit vectors e*^'\ whereas di'*^'-* = — r' a^paj^p^^ are 
the hopping integrals along the three directions a*^'^ of 
the next-nearest-neighbors, (t, t' are often called scal- 
ing parameter^iSl/) Following along the line of Ref.lSUfil^ 
one can calculate the variation of the band gap and of 
the energy at the K point under strain via hopping in- 
tegral modulation, and — after long yet not particularly 




FIG. 1; The two sublattices (circles and squares) of graphene 
and the three unit vectors e''' used in the text. </>, z are 
cylindrical coordinates of a tube, while ^ = -k /Q — 9c with 6c 
the chiral angle. 



insightful calculations — deduce ^ and in particular 



T3r' 



(7) 



The minus sign in the first equation is for electrons, the 
plus for holes. Note that there is no change in sign for v in 
going from holes to electrons. The other sign function is 
Si = — 1 (s2 — 1) where p is defined as a function of n, m 
as p= [{n — m) mod 3] . This difference in sign behavior 
between v and v' was previously recognized as causing 
indices dependence and non monotonicity in the doping 
induced shift of Raman frequencies and ano malous b ond 
contraction/expansion in zig-zag nanotube d^l^l^"*^^ . We 
will see that it controls chirality dependent expansions 
vs. contractions for hole or electron self-trapping. 

We assume that our fields only vary along the axial co- 
ordinate z, an ansatz corroborated by previous numerical 
calculations, as explained abov^^^. This treatment — we 
will see — is self consistent as it predicts very long po- 
larons. The hamiltonian density for the entire system is 
then 



n 



2m 



(8) 



where V^i is given by (|4|, E^p is given by (|6|, and c — 
27rr is the nanotube circumference. Minimization under 
normalization of ?/> with Lagrangia multiplier tp returns 
the system of equations of equilibrium for our system. 



4(/i+A)7„ = -i/'|V^|Vac 
4/47 - 4/3(g0S3 - qzCs.) = ~vc'i\i)'\^ loc 
4ilr]-\-4/3{q^C3 -\- qzSa) = -vssl'ipl'^ / ac 
4/3(703 -\- 77S3) -\- uj^qz = -v\i}\'^/eac 



4^(753 - '7C3) -u^q^ = 



'-diip +[u{q''/e + 7C3 -I- 7/S3) -I- v'-io] ip 



(9) 

which is linear in the elastic fields and can be solved easily 
for them in terms of |V'P- 

We find for the local deformation induced by an in- 
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jected electron or hole the equations 

-SpCsin(30c)|V'P 



(10) 



which show that strain is higher in regions where the elec- 
tron is localized. We find Co = j — . ^ , C = ^i^^^ 
inversely proportional to the circumference. The sign in 
the first equation is plus for electrons and minus for holes. 
The other two equations change sign depending on the 
chirality of the nanotube via Sp. 

From ^ we have for the total elongation and torsion 
of the tube AL = Jijo — 7) dz and rA0 = J rj dz, and 
thus from ( |10p and because of normalization of i/;, we 
obtain 



AL = + SpCcos(36'c) 
rAcj) = —Sp C sin(30c) 



(11) 



where, repetita juvant, the sign in the first equation is 
plus for electrons and minus for holes. Note that the 
overall elongation and torsion of the tube due to elec- 
tron/hole injection is independent of the actual shape 
of ip. Unfortunately, that also implies that self-trapped 
electronic excitations cannot be recognized by global ob- 
servations such has overall elongation, but solely via local 
mechanical features. 



Equations (111 generalize to arbitrary chirality what 



already found by Nisoli et al f or d oping of zig-zag nan- 
otubes^. The first result from (111 is that electron/hole 
injection induces no torsion for zig-zag nanotubes (which 
correspond to m = 0, thus = 0), at least in our approx- 
imation that neglects curvature corrections. On the other 
hand, torsion would be maximum for tubes approaching 
the armchair configuration (the armchair themselves are 
metallic and thus excluded from this study): that would 
be n = m + 1. Again from (111 we obtain that elongation 
is always positive for electrons in nanotubes with p — 2, 
whereas is always negative for holes in nanotubes with 
p = 1. The maximum in both cases is achieved by zig- 
zag nanotubes (9^ = 0). These results are in agreement 
with the density functional study of Ref.^**, who found 
elongation for electrons in (11,0) zig-zag, and contrac- 
tion for holes in (7,0) zig-zag. 

Analysis of the general case requires knowledge of the 
quantities ^o, If ^ > Co, the first equation of ( [TT| tells 
as that chirality can change the sign of AL Sp and thus 
predicts shortening for self trapped electrons, yet elon- 
gation for holes, in certain tubes. This seems to be the 
case, from a rough estimate: we use e = 1.42 A and the 
Harrison scaling for the hopping parameterJ^, which are 
r' = 2t', T = 2t, with t ~ 2.8 eV, t' ~ 0.68 eW^. For the 
speeds of sound in graphene we use vt = 1.4 10'* m s~*, 
wl = 2.16 10^ m s~* as in Mahan'^ ; finally from the den- 
sity of graphite 2.26 gm cm~'^ and the interlayer distance 
at 3.4 A the surface density of graphene can be estimated 
a ~ 7.6 10~^Kgm"2. We obtain ~ 2.5 10"^ A'^/2r, 
C ~ 8.7 10~^ A^/2r. For a larger nanotube of 1 nm in 



diameter, 10^ electrons can give a Inm elongation. Since 
the elongation and torsion are inversely proportional to 
the radius of the nanotube: these effects become stronger 
for small radii, for which our treatment is only a first or- 
der approximation. 

From (|9| we find that the inner displacement 



= 



-sp^{i + nm' 



(12) 



with /' = 4w^/a;p(e H- Z) ~ 0.35 A (computed from 
VT = 1.4 10'' ms-^ and wr = 3 10" s"*). Thus, the in- 
ner displacement is always parallel to the nanotube axis, 
while its absolute magnitude is independent of the chi- 
ral angle. On the other hand, its orientation depends 
on the chiral indices, thus suggesting that polarons can 
induce highly non monotonic anomalous bond contrac- 
tions/expansion and thus hardening/softening of optical 
modePaMHl. 

Finally, by substituting the expression for the elastic 
fields (10 12 1 into the last equation of (|9|, one finds a 



nonlinear Schroedinger equation for the envelope tp 



Cplp 



2m 



(32^ + 4x1^1 V) 



where the reciprocal length x is found to be 



4(A + A) 



(13) 



(14) 



where we have used the expression m — 2Trh^ / (Qect) for 
the effective mass of the electron obtained from the for- 
mula for the band gapf^. The well known self-bound so- 
lution of ( 13 1 



(15) 



and length 



corresponds to a polaron of energy €p = 
a = x~^- Remarkably unlike the amplitude of the elas- 
tic deformations, the length and thus the energy of the 
polaron does not depend on the chiral angle. Also po- 
laron length is the same for electrons and holes. This 
latter statement is in contrast with results reported by 
Verissimo-Alves^: while for the (11,0) nanotube they 
find about the same polaron length for electron (39 nm) 
and hole (40 nm), for the (7,0) case the electron po- 
laron has a length of 59 nm, while the hole is 21 nm 
long. We suspect that the discrepancy originates a sign 
error in the quantity in Ref citeVerissimoO defines as A^, 
as analysis of their Fig. 2 confirms. In our formalism, 
and using the Harrison scaling and parameters above, we 
have for a zigzag nanotube Xz — —3spt/2 =F 3t' and esti- 
mate Xz — —6.3 eV for electrons, -2.2 eV for holes, for 
a (10,0) nanotube-in the ballpark with the Ref.'^ val- 
ues of -8.3, -2.4 eV. Instead we find a sign mismatch for 
the (7,0) where our estimates are -t-2.2 eV for electrons. 
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+6.3 eV for holes, whereas they have -1.6 eV, +7.8 eV. 
The polaron length a = can be estimated from (14) 



to be of several hundred angstroms for nanotubes of a 
few nanometers in diameter, in agreement with estimate 
based on DFT resultM 

Finally, for completeness, we calculate the radial vari- 
ation of the nanotube Ar/r = p'^'^. From (10 1 and (15) 



we see that the maximum radial variation corresponds to 
the maximum of (at z = 0) or 



= ± 



k 

'2a 



2a 



cos(36'c), 



(16) 



positive sign for electrons, negative for holes. There is al- 
ways a radial shrinking in correspondence of self-trapped 
holes in nanotubes with p = 2, and radial expansion for 
electrons in nanotubes with p = 1 , as partially confirmed 
by density functional findings of Ref.^"*. 

While the casuistic of excitons is too rich to be dealt 
with her^i^l^ -^ve can still offer a few considerations for 
the Ell case. As couplings of electron and holes with 
anisotropic deformations 70 cancel each other, we must 



take 



and thus = in our equations. Hence, 



we can have optomechanical effects of both elongation or 
contraction, depending on the chiral indices, as predicted 
by pTl) with = 0. 



The treatment above applies to static polarons in not 



too small nanotubes, as we neglect the breaking of hon- 
eycomb symmetry brought up by the chiral vector. That 
would introduce new terms in the elastic energy and re- 
sult in curvature corrections for very small radiP^. Also, 
for small radii, orbital hybridization imposes a more so- 
phisticate tight binding treatment. For the propagating 
polaron, another symmetry breaking comes from the ve- 
locity vector, directed along the nanotube axis. Even for 
large nanotubes, we expects corrections (at lowest order 
quadratic) in the speed of the traveling polaron.!^. 

In conclusion, we have calculated the elastic deforma- 
tions and electronic structure of self-trapped electronic 
states in single wall carbon nanotubes. We found that 
elongation and torsion depends highly non monotonically 
on the wrapping indices, whereas the polaron length does 
not. The inner displacement is always axial and changes 
orientation depending on the chiral indices: this suggests 
a chirality dependent bond lengthening/contraction and 
optical mode frequency shift. Extension of this model to 
excitons is particularly interesting. 
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The expression for T4i is simphfied in the long polaron ap- 
proximation. Using the notation of Ref.l^, we can neglect 
terms in d^q-' (the dispersion in the optical branches) as 
long as {fj, — ii')/a^ , {\ — X')/a^ <S^ a, where a is the length 
of the polaron. Since {fi — ^'), (A — A') are the order of 
the square of the speed(s) of sound in graphen^^^or about 
10* m^s~^, whereas a can be expressed in terms of the 
graphite- like optical frequencjl^ a = cjp /4 ~ 10^® s~^ and 

10 nm, those 



we have for the length of the polaro: 
ratios are thus of the order of lO"** 



